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Abstract 

The algebraic and Hamiltonian structures of the multicomponent dispersionless Benney 
and Toda hierarchies are studied. This is achieved by using a modified set of variables 
for which there is a symmetry between the basic fields. This symmetry enables formulae 
normally given implicitly in terms of residues, such as conserved charges and fluxes, to 
be calculated explicitly. As a corollary of these results the equivalence of the Benney and 
Toda hierarchies is established. It is further shown that such quantities may be expressed 
in terms of generalized hyper geometric functions, the simplest example involving Leg- 
endre polynomials. These results are then extended to systems derived from a rational 
Lax function and a logarithmic function. Various reductions are also studied. 



1. Introduction 



One of the most studied integrable systems is the Toda lattice equation 

otot 

together with its various generalisations. Besides the inherent interest in such integrable 
systems, these Toda systems have a wide number of applications, form sorting theory to 
quantum field theory and differential geometry. The one dimensional Toda chain 

= ePn+i _ 2eP- + e^-i , (1) 
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may be written as a coupled system of difference equations 



P„,, = P(l-A-i)5„, 

where p„ = log P„ and A is the shift operator defined by A^p„ = pn±i , and the properties 
of such systems have been extensively studied, notably by Kuperschmidt 

This paper concerns a particular limit of these systems, the continuum or long wave 
limit. The idea is to look for solutions whose natural length scale is large compared with 
the lattice spacing, or, alternatively, one lets the lattice spacing tend to zero. In the limit 
the discrete label n becomes a continuous variable x and the operator A+^ — 2 + A~^ 
degenerates into a second derivative. Thus under this limit the two dimensional Toda 
lattice becomes 

dtdi dx^ ' ^ ' 

One intriguing property of these systems is that the notion of integrability is preserved in 
this limit, so (^) is in fact a multidimensional integrable system. In this paper this limit 
will be called the dispersionless limit, as mathematically it is the same as the limit which 
send the KdV equation into the Monge equation 



Ut = UUx ■ 

Many of the geometrical properties of these integrable systems have been derived in p|. 
One of the more recent motivations for the study of such integrable systems has come 
from their role in topological field theory [Q. 

Under this limit the difference system (0) becomes 

St = 
Pt = 

and in an earlier paper this system was 
to a set of modified variables defined by 



pi. w 

studied. This was achieved by transforming 



S = u + V , 
P = uv . 

In these new variables the Toda system @) becomes symmetric: 

Vt = VUx 

and various properties of this system, such as its associated hierarchy and Hamiltonian 
structure, were studied. In this paper these result are extended to the multicomponent 
Toda hierarchy. 



The multicomponent Toda hierarchy is defined in terms of a Lax function 



N~2 



c{p) = + p^s,{x, t) , t = {t,,t2,...} 



i=-l 



by the Lax equation 



(6) 



Here the bracket is defined by the formula 

' dp dx dx dp 

and, as usual, denotes the projection of the function O onto non- negative powers of 
p . More general Lax functions will be considered in section 6 and 7. 

Example 1.1 One obtains from the Lax equation with N = 2 ,n = 1 equation 
and, with = 3 , n = 1 the system (where C = p"^ + Sp + P + Qp^^ ) 

Pt = Qx, (7) 

Qt = ^QSx ■ 



The modified variables alluded to above are defined by the factorization of the Lax 
function 

1 ^ 

C = -l[[p + u,] (8) 

This transformation may be thought of as a dispersionless Miura map, and their use in 
the study of the dispersionless Toda hierarchy is due to Kuperschmidt Recall that 
the transformation from the KdV to the modified KdV equation comes from a similar 
factorization of the Lax operator for the KdV equation. These new variables Ui will be 
called the modified variables. 

There are a number of advantages in moving to these modified variables. Firstly, 
it puts all the fields on an egalitarian footing and any formula (such as those for the 
conserved charges and fiuxes and Hamiltonian structures) will be symmetric functions of 
these variables (explicitly, such formulae are invariant under the map Ui where 
a is an element of the permutation group). Secondly, it enables such functions to be 
calculated explicitly, using simple combinatorial arguments. In fact, rather than the 
calculations being of the form 

(1 + sum of N terms)" 

the calculations are of the form 



^(1 + single term) 

N 



and the symmetry between the modified variables means that only a single calculation 
has to be performed. The existence of such quantities and some of their properties follows 
from the general theory of these Lax equations (see, for example and the references 
therein). However, the explicit forms for these quantities are rarely given except, perhaps, 
for the lowest members of the hierarchy. The modified variables enables one to perform 
the general calculations with an arbitrary numbers of fields. 

One further property of these modified variables is best illustrated by means of an 
example. 

Example 1.2 In these modified variables equation becomes 

Ut = -u{-u^ + v.^ + w^) , 

Vt = ^v{+u.^ -v^ + w^) , (9) 
'>^t = ^w{+u^ + v^- w^) . 
where S = u + v + w,P = uv + vw + wu and Q = uvw . 

It is clear from these equations that one possible reduction of this system comes from 
the constraint v = w . In the original variables this constraint is 

4P^ + 27Q2 - 18PQS - P^S"^ + AQS^ = , 

this being the condition for the cubic equation C{p) = to have a double root. In 
the extreme case where all the fields are equal one obtains, after some rescalings, the 
dispersionless KdV hierarchy 

Ut„ = u'^Ux . 

Thus reductions are far easier to study using these modified variables. 



Such reductions are studied in section 5. 

This Toda hierarchy is related to the Benney moment equations which take the 
form 

dSn dSn+l , A OSq 

-7^ = ^ — + — , n = 0,l,.... 

ot ox ox 

The corresponding Lax function for a truncation of these equations at finite is given 
by 

^0 P 

and the corresponding Lax equation is 

^B,t = {('^b) + 7 ^b} ■ 



Example 1.3 For N = 2 the Benney system coincides with the dispersionless Toda 
equation (^). For N = 3 the Benney equations are (where Cb = p + S + Pp~^ + Qp^^ ) 

St = P'X , 

Pt = PS, + Q., (10) 
Qt = '^QSx ■ 

These equations possess an infinite number of conservation laws of the form 

g(n) ^ pin) 



one finds that 






- Ff^ 




- 

s 


Qf - 


fpin) 



and so the charge and fiux may be written in terms of a potential H^"^^ via 



dP ' dS ' 



the first few being 



jyw = p, 

^(1) = 2Q + 2PQ, 
fj{2) ^ SP'^ + QQS, 
^(3) = 12PQ + 12P^S + 12QS^ + APS\ 

With these one obtains equations satisfied by this potential, namely 

9P2 ~ 9^9g ' 

— ^ = p — ^ + 2g . . ■ 
9P2 9P(9g 

It is also possible to derive many other equations satisfied by this potential, for example, 

9g2 ~ 2~d§dP ~ 2~dPdQ~ ' 
Notice that if if*^"^ is a solution of the above equations, so is 

n + l dS 



This imphes that all lower conserved densities may be obtained from the higher ones by 
differentiation with respect to S. A similar observation for the KdV conserved quantities 
has been noted by Gerald Watts 0. There is an extremely simple formula for polynomial 
solutions for the potential H^'^\ It is simply given by the coefficient of in the expansion 
of 

p p^ 

i.e. as the integral 

2'Ki J 

These results generalize to arbitrary A^. 

In terms of the modified variables (where S = u + v + w etc.) this becomes 

Ut = u{vx + uix) , 

Vt = v{Wx + Ux) , (11) 

wt = w{u^ + v^) , 

and the equivalence between the Toda and Benney systems is given by the change of 
variables 

U = VW , V = wu , w = uv . 



With such a change the Benney system (|TT]) transforms into 



Vt = v{wu)^^, 

Wt = W{uv)a:, 

which is the N = 3 ,n = 2 flow of the Toda system (^. Without of use of modifled 
variables such an equivalence is less clear. 



Such changes of variable are discussed in greater depth in section 5. This will show the 
equivalence between the Toda and Benney hierarchies. An alternative way of showing 
this equivalence is given in section 6. 



Throughout this paper the binomial coefficients y ^ ) should be interpreted in terms 
of F-functions, i.e. 



Via 



r(a-6 + i)r(6 + 1) ■ 



This enables the binomial coefficients to be defined for fractional values of a and b . Such 
coefficients may also be defined for negative values if one manipulates the formula formally, 
using the basic definition r{z + 1) = zr{z) . For example 

f-2\ ^ r(-i) 
v+2; r(-3)r(3)' 

(-2)r(-2) 

r(-3)r(3) ' 

(-2)(-3)r(-3) 

(+2)(+i)r(-3) • 

For notational simplicity binomial coefficients have been used, rather than the more cor- 
rect F-function notation. 



2. The multicomponent Toda hierarchy 

It follows from the general theory behind such Lax equations [|l], ^ that the quantities 

QW = -L/£i^^ (12) 

ZTTZ J p 

are conserved with respect to the evolutions defined by the Lax equation (P). 
Proposition 2.1 

The conserved charges defined by ([T^) are given by the formula 

Proof 

It follows from (||) that 



i=l 



i=i y 



and so, using the binomial expansion, 

N oo 



= n E V <'p 



N-1 \ ^,^:-^r,-''"i 

i=l ri=0 

Using the residue theorem the integral formula (p!^ ) may be easily evaluated, 

N oa . n \ 

Q(") = coefficient of p"" in H E P^'M ' 

{ ^ 

= E n 



i=l 



N-i ) 



Hence the result. 



□ 



Example 2.2 Let = 2 . Then 



Q 



in) 



E 

•+s 
n 

E 

r=0 



n\ n 



r+s=n 



n / X 2 

n 



r \S 



U V 



This was calculated, using a different approach in 0]. 



These charges are conserved with respect to all the times, i.e. 

dtn ~ dx ^ ' 

for some function A*^"^'"-' . Before this function is calculated it is first necessary to find the 
explicit form of the evolution equations for the fields Ui . 

Theorem 2.3 

The Lax equation (P) implies the following evolution equations for the fields 



Ui 



(n) 



ij ^j,^ 



where 



TV 

n 

k=l 



n 



N-1 



Ti - 1 



and 

B 



[n) _ n MT^ 

~ N -1 ^ 



AT 

n 

k=l 



N-1 ] 



Tk 



n 



N-l 



U, 



r 



3 



(15) 



' . (16) 



Proof 



Since 



N 



we have, on differentiating. 



N 



dt 



N 



P j=0 I i=l 

N 



U 



j,tn ' 



Ui 



"lit 71 



and similarly 



dC 

dx 

dC 
dp 



N 



t^lP + Ui 



1 A 1 



P t^lP + Ui 



The first of these derivatives implies that 

= Res [C-'{M,C}], 

P=-Ui 

where, for convenience, we have defined M. to be the quantity M. — {jC^^)+ . The explicit 
form of M. is easy to calculate. 



So 



N oo , _n_ 

cw^i = n E ' 



1=1 ri=0 

oo 



Ui'P \ 



E 

t=0 r .V^iV 



{'•i •Ei=i '•i=*} 



TV 

n' ^-^ 



i+ 5 



p^E E 

t=0 r .V^JV 



iV . n 

n ^-^ I 



The derivatives of Ai are more complicated, but are easily evaluated. Using these we 
obtain 



n 

Ui^t = Res Ep' 



^ / -0—\ 1 r ^ 1 r ^ 1 

n(^-i)<^ y(n-t)ui,^—+[i-y^—\ 



and hence 



u. 



iitn 



(n) 



where 



A 



in) 



i=0 



(n) 



E(-i)""' E 



k= 

N 

n 

k=l 



Tk 



n \ / n 



To obtain the explicit forms (|T5D and (|T6D of these one has to perform a comphcated 
resummation. In the simplest (A^ = 2) case this amounts to using the identity 



E(-i)"^' E 



n\ n 



r J \s 



n+r~t 1 
SU V 



n 



E 



t=Q r+s=t ^ ' ^ ' 1 

the proof of which has used various binomial identities. 



n — 1\ ( n — \ 



r s 
U V 



Example 2.4 For general N the first {n = 1) flow is 



N 



iV - 1 ^ 



Ui 



and the (A^ - I)**' flow is 



Ui,t = Ui[ Y[u 
For N = 2 these are the same. 

Example 2.5 With = 3 and fields u , v and w the above formulae give: 

Ut = Aux + uBvx + uCwx 

and cyclically, where 



□ 



Br, 



Cn 



,n 



--1) y , 

p+q+'r=n ^ ^ / \ i 



- - 1 

2 



n \ / 71 



2 W 2 1 ^,P^;9y;'- 



2 ^ 

p+g+r=n— 1 



p 



2 ) , 



E 

p+g+r=7i— 1 



- - 1 

2 



P / \q 



r 



Example 2.6 A similar calculation may be performed for the Benney hierarchy. In terms 
of the modified variables the first fiow is 



this being the multicomponent generahsation of (]TT|). 



Having found the exphcit form of the evolution equations one can now derive the flux 
/\ ("!.>") corresponding to the conserved charge (P^P^ . 



Proposition 2.7 



/V - 1 ^ 



m + n 



i=l jj^i 



where 



p{m) _ ^(m,l) 



E 



{'■fe •Z]fc=i '■'=='"+1} 



N-l \ yj'k 



Proof 

It follows from the conservation law (|l^) that 



N 



in). 



N 



'Ui,t 

1=1 1=1 



EA (n,m) 



Mi 



and using the evolution equation this becomes, on equating the various coefficients of the 
derivatives, 



JV 

I 

1=1 



JV 



i=l j^i 



It is clear, by power counting, that A*^™'"^ must be a homogeneous function of degree 
(m + n) in the fields, so by Euler's theorem 



N 



1=1 



Since the Ai and -Bjj are known it is now just a matter of substituting these into the above 
formula and simplifying the results. In fact 



E ^ 



N 

n 

A;=l 



^-1 ul" 



N-l 



1 



(17) 



the proof of which, once again, involves the use of various binomial identities. Let F^^"* = 
^(m,i) ^ ^Y^e flux corresponding to the first time fiow. The above formulae simplify to 
give 

r / ™ \ 



so 



1 



A^- 1 
1 



(m) n(m) 



AT 

n 

fc=i 



N 



by using equations (|r^). Thus 



AT - 1 ^ 



m + n 



The first order flow for general N has been calculated in an example and it follows from 
this that 



Inverting this gives 



and so 



u 



.r){n) _ p(n) 



/V - 1 ^ 

where Gj,- is the matrix with zero entries along the diagonal and — 1 everywhere else. 



□ 



Note that A^™'") = A^'"-'") , i.e. 



It follows from this, together with the commutativity of the flows, that 

d d ^ 



^ {m,n) 



df df 



for some function F . This is called the free energy and is of great importance in topological 
field theory ^. There exist various algebraic relations between these A^"^'") , these 
being derived from the differential Fay identities 0, ||, |^. The matrix Gij is very closely 
connected to the components of the metric g which defined the Hamiltonian structure of 
the hierarchy (see section 4). Where this is accidental or indicative of some deeper result 
is unknown. 



3. The algebraic structure of the charges 



Consider the conserved charges for the N = 2 hierarchy 



Q 



in) 



E 

r+s=n 



n\ n 



If one singles out one of the fields, say u, this may be written as 



Q 



(n) 



r=0 



and so the properties of these charges can be reduced to the study of the polynomials 



r=0 



These polynomials are examples of hypergeometric functions, 

fn{x) = 2Fi{-n, -n, l;x) 

(since n is an integer this series will automatically truncate leaving a polynomial of degree 
n), and so the conserved charges may be written as 



Q 



(n) n 



2Fi( - n, -n, 1; 



Using the identities 



Pn{z) = 2Fi{n + l,-n,l- 



1 - z 



(where Pn{z) is the Legendre polynomial of degree n) one obtains 



Q 



iu - vrp„ 



U + V 



u — V 



and so a generating function for these charges may be constructed using the well known 
generating function for Legendre polynomials. Such a generating function was obtained 
in [Q using a different approach based on the construction of recursion relations. 

In the multicomponent case one may perform a similar calculation and express the 
charges in terms of generalized hypergeometric functions |jlO[- Defining new variables Xi 

by 

{ UN, t= N , 

i ^ N 



Xi 



UN 



equation ([T3|) may be written as 



n 



N-1 

n \N 



N 



E 



X 



r,z--n^=ir(i + Mr)+7fc)' 

(11 \ 
^^^j xl F(/i,7;x) 



where, for convenience, is defined to be 



r n ^JV— 1 

— ... X 7v_]^ , 



and Ti are the components of the vector r . The hnear map fi : Z 



N-l 



T?^ is defined by 



iV-l 



■3' 3 



i=i 



where 



+1 


.. 


• ° \ 





+1 .. 


. 





.. 


. +1 


+1 


+1 .. 


. +1 


-1 


.. 


. 





-1 . . 


. 





.. 


. -1 


-1 


-1 . . 


. -1/ 



and the constant vector 7 is 



7 








n 



\N-l) 



N-l 

n 

N-l 



\ n J 

These manipulations are a manifestation of the Ore-Sato Theorem for generahzed hy- 
pergeometric functions. Thus the charges are defined in terms of a hnear map and a 
constant vector 7 , or in terms of a lattice B inside 7?^ as the image of Z^~^ under the 
linear map /i . These generalized hypergeometric functions a related (but not identical) to 
those studied by Gelfand et al. ||10[ . 



The fluxes (|T8D also have a similar description. 



Vat- 1^ 



n N 



xl F(/i,7';x) 



where the linear map /i is the same but the constant vector has changed to 





7 =7 + 



1 





So both the charges and the flux are defined in terms of the same lattice B G 7j , the 
difference being in the constant vector. One remaining problem is to understand the 
structure of the more general fluxes A*^™''") in terms of these structures. 



4. The Hamiltonian structure of the hierarchy 

The evolution equations studied in this paper are speciflc examples of hydrodynamic 
equations, i.e. they are of the form 

«i = E^^("K' (19) 

the simplest non-trivial example being the Monge equation Ut = uu^ ■ One must, in 
general, take care with the indices on the flelds as they transform as tensors under an 
arbitrary transformation Ui 'Ui(u) . However, to avoid expression such as [u^Y^ all 
indices have been lowered for notational convenience. There is an extensive literature on 
such hydrodynamic equations (recently they have attracted much interest from their role 
in topological fleld theories). One important result concerns the Hamiltonian structure of 
these systems. 



A system ([T9|) of hydrodynamic type is said to be Hamiltonian if there exists a Hamil- 
tonian H = J dxh{u) and a Hamiltonian operator 

A^^ =gHn)^ + b\{u)u': (20) 
which deflnes a skew-symmetric Poisson bracket on functionals 



which satisfies the Jacobi identity and which generates the system 

ui = {u\x),H}. 

Dubrovin and Novikov proved necessary and sufficient conditions for A''^ to be a 
Hamiltonian operator in the case when g'^^ is not degenerate. These are: 

a) g = {g^-')^^ defines a Riemannian metric, 

b) V'\ = —g^^V\^ , where F-'^;, is the Christoffel symbol generated by g , 

c) the Riemann curvature tensor of g vanishes. 

The system (|T^) may then be written as 

«i = {g''V sV ,h)v^. (21) 

where V is the covariant derivative generated by g . Expanding these in terms of the 
Christoffel symbols gives 



where 



+ X! u^BikU^ (no sum) 



A' = Y.9'^didjh-Y.g'^Ti^d,h, (22) 

3 s J 

B.k = Y.9''dkd,h-J2g^^Ti,d,h (k^z), (23) 



s,3 



Thus to find the Hamiltonian structure for the hierarchy (^) one needs to find the metric 
and the Hamiltonians. 

One approach is to first diagonahse the system, there being the following simple for- 
mulae for the metric coefficients for diagonal systems |[T2| : 

- -^^\og g^j 



where VJ{u) = l^(u)(5* (no sum). For N = 2 such an approach was adopted in [Q. 
However, to find the corresponding transformation for arbitrary N is difficult and another 
approach is needed, the N = 2 result suggesting the form of the general result. 

Theorem 4.1 

The Hamiltonian structure for the hierarchy (|^) is given by (pO]) , where the zero- 
curvature metric is given by 

and the Hamiltonians generating these fiows are given by 



n 



(n) 



where the Q^"^ are conserved charges given by (|T2|,|T3[) . 
Proof 

This g defines a non-degenerate Riemannian metric, and by changing coordinates to 
Ui = log(Mi) it is clear that the metric has vanishing Riemannian curvature tensor. The 
components of the inverse metric are 



-{N-2) UiUi 

UiUj 



(no sum) , 



The only non-zero Christoffel symbols are F^^ (no sum). This is easily proved by noting 
that in the Ui coordinates the Christoffel symbols are all zero. The Christoffel symbols 
can then be calculated using transformation law for these symbols. This yields 

V\i = -— (no sum) 

Ui 

with all other components being zero. Thus, by the theorem of Dubrovin and Novikov, 
this metric defines a Hamiltonian structure. To show that this is the dispersionless Toda 
hierarchy one has to verify that the Ai given by (^2|) and the Bij given by (^3|) and 
/i*^") = ^Q*-"-* are the same as given by equations (0) and (|16D in theorem 2.3. 



Ui 



Now 



n ^ 



and substitution of this into the above gives 



N 

n 

■t=i 



N-i \ yj;^ 



A 



(n) 



n 



N -I 



N 

n 

L k=l 



N-1 \ 



7V-1 
r,; - 1 



Ui 



The proof that b\^^ = is similar, both using various binomial identities. Hence the 
result. 



□ 



This theorem thus determines the Hamiltonian structure of the dispersionless Toda 
hierarchy. 



5. Miscellany 



This section contains some miscellaneous results on the multicomponent Toda hierarchy. 



5.1. Reductions 



Besides the discrete symmetry Ui i-^ 'Uo-(j) , where a G Sn is an element of the permutation 
group, the muhicomponent Toda hierarchy possesses many other symmetries, as befits 
an integrable system. For example, one has the continuous symmetries corresponding to 
time and space translations, and various scaling symmetries. These may be used to derive 
examples of integrable dynamical systems. 



With the ansatz 



Ui{x,t) 



(nx 



for the n**^ flow (this being a particular example of a scaling symmetry) the x-dependence 
cancels and one is left with a dynamical system for the Zi variables. 



Example 5.1 With n 
system 



1 the above ansatz reduces the system ([TT|) to the dynamical 



Zl{Z2 + Z3) , 

22(^3 + Zl) , 

ZsUi + Z2) . 



This system has been studied and integrated by Bureau [O . 



In general, such dynamical system are of the form 

Zi = Zi (polynomial in the Zj variables) 

(since A]^ given by equation (|15]) always has Hi clS cl factor), which suggest on should 
introduce variables 0i = log Zi . The resultant dynamical systems are then in the so-called 
Lotka-Volterra form, the nonlinear terms being sums of products of exponentials. 

As mentioned in section 2, one may obtain other systems by a conflation of the fields 
(or equivalently, by having multiple roots in the equation C{j)) = 0). As an example, the 
system (^) 



1 / N 

Ut = -u[-u^ + v^ + w^) 
wt = -w{+u^ + - w^) 



becomes, one setting v = w . 



Ut = -u{2vx - Mx) , 

1 



{1,1,1,1} 



{1,1,1,1,1} 



{2,1,1} 



{2,1,1,1} 



{3,1} {2,2} 



{2,2,1} {3,1,1} 





{4,1} {3,2} 




Figure 1: Possible reductions of the = 4 and N = 5 Toda hierarchies 



and finally, on setting u = v , 



Ut = -UUcc . 



The Hamiltonian structure of these later systems is obtained by restricting the metric 
(which defined the Hamiltonian structure of the original system) in the first case to the 
surface v = w and in the second case to the line u = v = w . Schematically these systems 
may be written as {1, 1, 1} , {2, 1} and {3} respectively to show the confiuence of the fields 
(the discrete permutation symmetry ensuring that {2, 1} is the same as {1, 2}), and hence 
the possible reductions may be summarised in the diagram 

{1,1,1} ^{2,1} ^{3}. 

With more fields the are many more possible reductions and a simple combinatorial 
argument gives the number of such subsystems (including the original system) to be 
2iV2 + 7+(-l)^ 

for the A^-component Toda hierarchy. Figure 1 shows the possible 

8 

reductions for = 4 and = 5 . 



5.2. Special flows 



It is clear from the Lax equation (^) and the integral for the conserved charges (|T2|) that 
if n = (A^ — l)m for some integer m then some simplification might occur, since for such 
a value of n the Lax equation 



dC 

dim. 



{(£-)+,£} 



(where im = t(N-i)m ) and the integral form for the conserved charges 



2m 



P 



do not involve fractional powers of the Lax function. These flows are best studied by a 
change of variables 



Ui 



n 



this being the generalisation of the one used in example 1.3. In these variables the 
evolution equations for these flows take the form 



u 



(25) 



where 



A. 



(m) 



B. 



(m) 



(m 



m 



1) E 



N 



( m 



m — 2 

n-l 



E 



N 

n 

L k=l 



and the corresponding conserved charges are 



m 

Tk 



u 



m — 1 



Ui 



m — 1 



u 



'j ■ 



Q 



(m) 



E 



N 

n 



m 



Ui 



Superficially these equations look very similar to their unhatted counterparts. However, 
it is important to note that only appears in a subordinate role, i.e. instead of the 
summation being over the set {rj = (^ ~ l)m} the change of variable enables this 

to be written as a summation over the set {r^ : J^f^i = "^j irrespective of the value of N . 
Similarly does not enter the binomial coefficients that appear in the above formulae. 
One consequence of this is, if one denotes the evolution equation (p5D by r(A^, m) , then a 
possible reduction is to set un = and the resultant system is just the (A^ — l)-component 
system, i.e. 



T{N,m] 



for all values of m . 



The Hamiltonian structure of this system is easily derived by performing the change 
of variable on the metric (p^), the result being 



{N-iy 



1=1 



Ui 



Uj^ Uj^ 



(actually, the inverse metric 

1 otherwise, 

is simplier). 

These calculations are not, perhaps, very illuminating, and the motivation for this 
subsection is best seen by means of the following examples. 

Example 5.2 For general N the (A^ — 1)*^ flow (i.e. m = 1) has been calculated in 
example 2.4, namely, 

Ui^t = UiiYluA . 

In the hatted variables these simplify to 

This is the first flow of the Benney hierarchy (see example 2.6) written in terms of the 
modified variables. Note that this contains only quadratic nonlinearities irrespective of 
the value of , unlike its unhatted counterpart whose nonlinearities depend on the value 
oiN . 

Example 5.3 For = 3 the Lax equation (|^) gives the evolution equations for n = 2 

Ut^ = uwVx + uvWx (and cyclically) 

and for = 4 

Ut^ = 2uvw{v + w)Ux + 2uw{2uv + uw + vw)vx + 2uv{uv + 2uw + vw)wx (and cychcally) 
together with the conserved charges 

Q*-^-* = uv + vw + wu , 

_ _^ ( a\ 00 OO 00 / \ 

' = uv+vw+wu + Auvw{u + V + w) . 
In terms of the new hatted variables u = vw etc. these flows become 

= u{v^ + w^) , 

= {2uv + 2uw)ux + {2v? + 2uv + Auw)vx + {2u^ + Auv + 2uw)wx 
(together with their cyclic permutations) and the conserved charges become 

Q^"^^ = -u^ + -0^ + + 4:{uv + vw + wu) . 

These flows are thus simpler in the hatted variables. Note that it is now possible to set, 
for example, w = and these formulae reduce to those for N = 2 (for which there is no 
distinction between hatted and unhatted variables). 



6. Rational Lax functions 



The Lax function so far studied in this paper has been, apart from an overall p'^ , a 
polynomial function. However, the Lax formalism is valid for much more general functions 
such as rational functions: 



{P + Ui 



T-fN+M 
lli=iV+l 



(26) 



or 



M+N 

>c= n ip+^iY' 
1=1 

provided N > M . Here = ±1 but this can be extended to cope with repeated roots 
and zeros. The N > M condition enables this rational Lax function to be written as 

N+M 

C = polynomial of degree {N — M) + ^ simple poles . 

i=N+l 

It should be clear from the earlier proofs that they may be extended to such rational 
potentials with little change. The results in this section will therefore be given without 
proof. 



The conserved charges, now defined by the integral formula 

27ri J v 



are 



N+M , e^n 



n 

i=l 



N-M 



and the evolution of the fields, given by the Lax equation 



are 



where 



A 



(n) 



N+M , e^n 

n ~ 

L k=\ 



N-M \ yj'k 



N-M 



and 



N-M 



N+M , efcn 



n 

fc=i 



N-M \ y^'h 



N-M 



N-M 



Example 6.1 The simplest case is A*" = 2 , M = 1 , and the Lax function will be written 
as 

^ ^ {p + u){p + v) 
p + w 

Hence a consistent reduction is to set w — , with which the results simplify further. The 
conserved charges are 

The binomial coefficients have to be interpreted formally, as indicated in section 1. For 
example, the first two charges are 



Q^^^ — u + V — w , 

Q^"^^ — + + 3w'^ + 4:{uv — vw — uw) , 



and the first two flows are given by 




/O 







(u 


:: 




V 


-2w J 




I w 



and 



2u{v — w) 2u{u + V — 2w) 

2v{u + v — 2w) 2v{u — w) 



—2u{u + 2^ — 3tu) 
-2v{2u + i; - 3w) 




2w{u -\- 2v — Zv) 2w{2u-\-2v — Zw) —Qw{u-\-v — 2w) 

Once again, care has to be taken in evaluating the binomial coefficients. 

Example 6.2 Let M = iV - 1 and = for i = + 1 , . . . , + M . With these the 
Lax function takes the form 



Thus the Lax function Cb for the Benney hierarchy is a special case of the rational Lax 
function. Thus one obtains an explicit family of conservation laws and evolution equations 
for the Benney hierarchy. 



In fact, the rational Lax function (^) itself may be considered as a reduction of the 
infinite component Lax function 

N'-l 

Coo E S^P''. 

n=—OD 

By formally expanding (|2^) one obtains the above Lax function, with various conditions 
of the fields. This and some other reductions of this Lax function, and the resulting 
systems, have been studied by [Q. 



7. Logarithmic charges 



In section 4 it was shown that the conserved charges (for N = 2) may be written in terms 
of Legendre polynomials 

QH = («-.rP„(^). (27) 
This suggests that one should introduce the new variables 



r 

cos^ 



u — V , 

U + V 

u — V 



In these new variables the differential equation for the conserved charges 

{uQu)u = {vQv)v 

becomes the axially symmetric Laplace equation, and so one obtains the conserved charges 
( P7| ) . However there is another family of charges given by Legendre functions of the second 
kind, the first few conservation laws begin: 



- log{uv) 
''-(« + v) \og{uv) 



-{u + v) 
-uv log{uv) + - {u'^ + f ^) 



This raises the question of how these charges are related to the Lax formalism. The 
answer is that they are given by the formula 



Q^-'^ = ^{C- {log C-c^)^ 
2m J p 



(2J 



where 



n ^ 

Cn = E ~ ' Co = . 

j=l ^ 



This results in the scaling symmetry |]T5|, |T6 

d 



dC 



[/:"(log/: - c„)] = n[C^~\\ogC - Cn-i 



since Cn = c„_i + 1/n. Once again using the modified variables this integral may be 
performed for all values of n and . For simplicity we restrict attention to the N = 2 
case. Conceptually the general case is the same, just notationally more complex. 

Theorem 7.1 

The logarithmic charges defined by (ESI) are 



Q(") = I \ogiuv) E ( ; ) ( ^ ) u^v' - E ( ^ ) ( ^ ) s)u^v^ (29) 

r+s=n ^ / ^ / r+2=n ^ ^ ^ ''^ 

where 

r r+s 

j=l J j=s+l J 

and c(0, r) = (so c(r, s) = c(s, r) for all r and s). 
Proof 

Care has to be taken with the interpretation of log£ to avoid logp-terms. Since 

log £ = + logp + log ( 1 + - ) + log ( 1 + - ) 

\ p/ ^ 

and 

log£ = - logp + log f 1 + -) + log f 1 + -) + \ogiuv) 

then 



log C = - 

^ 2 



^ log (l + -) + log (l + -) + \og{uv) + log f 1 + -) + log f 1 + - 



p/ ^ p 



To evaluate the integral one Taylor expands the function as power series in p and uses 
the residue theorem. Since 



jr-iog(i + -)=p^yyy^r) r) L^u^+Wys -(r+s+w) 



2vr^J p t+1 

The other parts of the integral may be evaluated in a similar way and so the final result 



IS 



n \ I n 
r J \ s 



r+s=n 



Note that there is an ambiguity in the definition of Q^"-* , since Q^*^-* + XQ^^^ is also 
a conserved charge. The numbers c„ have the property that the coefficient of (m" + f ") 



in the terms not involving log(-uf ) in the above expression is zero and this removes this 
ambiguity. Exphcitly this coefficient is 



E 

r+t=n— 1 



n 



r J t+l 



and this is, using the identity 



1)* / n \ " 1 



J^\m + t) t+l \m- l) ^ j 



n > m . 



zero, as required. Rearranging and resuming (using the above identity) results in the 
expression (pH]). 



□ 

These logarithmic charges define a new set of evolution equations defined by the Lax 
equation 

dC 

— = {[r{\ogC-Cn)U,C}, n = 0,l,2,... , 
or, equivalently, by the Hamiltonian equations 

with if*^") proportional to the logarithmic charges Q*-"-* . These fiows all commute, both 
with themselves and those defined by (^, and the corresponding Hamiltonians are in 
involution with respect to the Hamiltonian structure defined by (PP]). The techniques 
developed in sections 2 and 4 may easily be extended to derive these more complicated 
evolution equations explicitly. However, for = 2 one may sidestep these calculations. 
It is straightforward to show that any fiow that commutes with the basic equation 
must be of the form 



Ur = FUx + uQVx , 
Vr = FVu + VQU^ , 

where F and Q are any functions of u and v (not necessarily symmetric) which satisfy 
the equations 

Fu "vQv 1 (30) 

or equivalently, the conservation law Qt^ = F^ . Thus to find the evolution equations 
associated with the logarithmic charges given by (p8]) one only has to solve the above 
equations for the fiux F . 



Lemma 7.2 The logarithmic charges (p8|) give rise to the evolution equations 



= F'-^^v^ + vQ^'^^u,, 



where the Q^"^ are the charges defined by ( P^ and the F*^") are the corresponding fiuxes 
with respect to the ti fiow, Exphcitly 



1 



= - \og{uv, _ 

^ r+s=n+l V ' / V 



E 



n\ n 



c(r, s)u^v^ 

r+s=n+l 



where 



(31) 



n \ / n 
r J \ s 



c r, s 



1 + 



1 + 1 \ + 1 

2(1 + n) [ r )[ s 



Proof Straightforward. One just has to integrate equation 



Example 7.3 The simplest such logarithmic fiow is 



■"TO 



or, in the original variables, 



^(m + V)U^ + \l log(Mf ) , 
1 1 

- (m + f )f a; + -V log(Mt;) , 



□ 



S. 



T"0 



TO 



^S-^. + ^logPF,, 
^^P, + ^logP^,. 



Calculations of this type were first performed (for N = 2 and small values of n) in 
15], p!6| and then (for = 3 and small values of n) in [O. 



8. Comments 



The results of this paper stem from the simple observation that by factorizing the Lax 
function the calculations simplify dramatically. As an example, this enables the metric 
which defines the Hamiltonian structure to be found for arbitrary . Even for A^ = 3 
the metric, when written in terms of the original, unmodified, variables is somewhat 
unmanageable and the degree of complexity increases rapidly as A^ increases. Not all of 
the possible formulae have been calculated in this paper, though it should be clear how 
these may be derived. For example, one may regard the r„ fiows as negative tn fiows (i.e. 
t-n = Tn), and so one may derive explicit forms for the function A*^"^'") for all m , n G Z . 
Similarly, the Hamiltonian structure for the evolution equations derived from the rational 
Lax function (^) has not been calculated; to do this should be straightforward. 



The use of generalized hypergeometric functions to study the conservation laws follows 
from an earlier result where it was shown that the N = 2 conservation laws were 
expressible in terms of Legendre polynomials. This results was used to construct new 
families of integrable (or at least soluble) hierarchies with similar properties to the original 
Toda system, namely 



where a and c are free parameters. The conservation laws of this system depend on the 
hypergeometric function 2Fi{—n, n + a — l,c; z) . One interesting possibility is to use the 
properties of the generalized hypergeometric series to perform a similar calculation for 
higher values of N . 

As remarked at the very beginning of this paper, these systems come from the dis- 
persionless limit of integrable dispersive systems. The Lax formalism for such dispersive 
systems comes from the replacements 



p ^ A, 
{f,9} — ' [f,9], 

where A is the shift operator introduced in section 1 and [/, g] = fg — gf where now / 
and g are operators depending on A . Alternatively, one may use a deformation of the 



Poisson bracket [|18|, this having the advantage that one retains a geometrical description 
of these hierarchies . While all this clearly works for polynomial Lax functions, there is 
a problem of how to find dispersionless analogues of the systems derived from the rational 
Lax function (p6|), the difficulty coming form the interpretation of the operator (A + -u)~^ . 
One may to resolve this problem has been recently proposed in . 



Finally, it must be stressed again that the existence of the results in this paper have 
been known for many years. However, apart from some trivial examples, the precise forms 
of these results have not been calculated. We have adopted an utilitarian view that such 
implicit formulae should be calculated explicitly, and the use of modified variables, as 
advocated in this paper, provides a simple way to do this. 
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